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We present a comprehensive calculation of the energies, splittings and electromagnetic decays 
of the low-lying levels of the bottom charmed meson system. In order to incorporate running 
coupling constant effects, we choose Richardson's potential for the central potential and take the 
"q^ ' spin-dependent potentials from the radiative one-loop calculation of Pantaleone, Tye and Ng. The 

, effects of a nonperturbative spin-orbit potential are also included. Our parameters are determined 

■ from the low-lying levels of charmonium (avg. dev. of 19.9 MeV) and of the upsilon system (avg. 
dev. of 4.3 MeV). We carry out a detailed comparison with the calculations of Eichten and Quigg 
and the lattice calculations of the NRQCD collaboration. Our predicted result for the ground 
state energy is 6286lg 5 MeV. Our results are generally in agreement with the earlier calculations. 
However, we find the two lowest 1 + states to be very close to the j-j limit, in contrast to some of 
the earlier calculations. The implications of this finding for the photon spectra of the IP and 2S 
states are discussed in detail. Some strategies for the observation of these states are discussed, and 

ff) \ a table of their cascades to the ground state is presented. Our calculated value for the ground state 

lifetime is 0.38 ± 0.03 ps, in good agreement with the recent CDF measurement. 
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■ The CDF collaboration M has reported the discovery of the B c system in 1.8 TeV p — p collisions at Fermilab. 
They have observed about 20 decays in the J/i/>-lepton channel, which are interpreted as decays of the ground state. 
For the mass of the ground state, the CDF collaboration quotes Mb c = 6.40 ± 0.39 ± 0.13 GeV. Their value for the 
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I. INTRODUCTION 



lifetime is tb = 0.46+q }es i 0-03 ps. This state should be one of a number of states lying below the threshold for 
Oh! emission of B and D mesons. Because these states can not decay by gluon annihilation, they should be very stable in 
comparison with their counterparts in charmonium and the upsilon system. 

The purpose of this report is to give a detailed account of the energies, splittings and electromagnetic decay rates for 
the B c states below the continuum threshold and to propose strategies for detecting some of these states. We will use 
a potential model that includes running coupling constant effects in both the central potential and the spin-dependent 
potentials to give a simultaneous account of the properties of charmonium, the upsilon system and the B c system. 
We choose Richardson's potential to represent the central potential and insist upon strict flavor-independence of 
its parameters. Since one would expect the average values of the momentum transfer in the various quark-antiquark 
states to be different, some variation in the values of the strong coupling constant and the renormalization scale in 
the spin-dependent potentials should be expected. In order to minimize the role of flavor-dependence, we use the 
same values for the coupling constant and the renormalization scale for each of the levels in a given system and 
require that these values be consistent with a universal QCD scale. Since one can calculate the central potential and 
the spin-dependent potentials from first principles on the lattice ||, it should be possible at some point in future 
calculations to connect subsets of these potential parameters in more fundamental ways. Such a program would be 
an important step towards finding a more rigorous way to formulate a potential model calculation. 

In 1991 Kwong and Rosner [Q predicted the masses of the lowest vector and pseudoscalar states of the B s and 
the B c systems using an empirical mass formula and a logarithmic potential. Eichten and Quigg ^ gave a more 
comprehensive account of the energies and decays of the of the B c system that was based on the QCD-motivated 
potential of Buchmiiller and Tye [|| . Gershtein et al. Q| also published a detailed account of the energies and decays 
of the B c system and established contact with QCD sum-rule calculations. Both of these latter calculations included 
running coupling constant effects in the central potential, but used spin-dependent potentials that were restricted to 
the tree level. 

One of the most important goals of the present calculation is to extend the treatment of the spin-dependent potentials 
to the full radiative one-loop level and thus include effects of the running coupling constant in these potentials. In 
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a previous calculation ]|| we have shown that such effects offer a substantial improvement in the calculation of the 
spectra of charmonium and the upsilon system. In particular, it becomes possible to offer a good account of both 
the fine structure splittings and the hyperfine structure. Our previous calculation used the renormalization scheme 
developed by Gupta and Radford [p|JIo| . However for this calculation, we have chosen the modified-minimal subtraction 
scheme used by Pantaleone, Tye and Ng (PTN) to extend their perturbative QCD calculation to the one-loop 
level. We supplement the PTN formalism with a long-range spin-orbit potential to be consistent with the Gromes 
consistency condition |l2j and recent lattice calculations || . The flavor dependence allowed the string constant in the 
1991 calculation is not permitted in the present calculation. Thus this calculation should also be viewed as part of 
the effort to formulate potential model calculations under a more restrictive set of assumptions. We note in passing 
a number of additional important calculations of the properties of the B c system |l3|-|l7| . 

In Sect. |n| we present preliminary calculations of the lowest vector and pseudoscalar states of the B c system with 
several potential models. These exercises are similar to some carried out by Eichten and Quigg Q and represent one 



way of estimating the uncertaintities in the predicted energies of the new system. In Sect. Ill we introduce Richardson's 
potential and the expressions necessary to include the one-loop corrections in the spin-dependent potentials and discuss 
the determination of our parameters from charmonium and the upsilon system. The formalism required for the mixing 



of the P states in the B c system is presented in Sect. IV. The last section contains our results for the energies and 



decays of the low-lying B c states and a discussion of their implications for the development of strategies to detect 
these states. 

II. SOME POTENTIAL MODELS 

The three systems that we consider in this paper, upsilon, charmonium and B c , are often considered as nonrelativistic 
systems, and thus our treatment is based upon the Schrodinger equation with a Hamiltonian of the form, 

H = p 2 /2 f i + V(r) + V SD , (1) 

where [i = m\m,2i '{mi + m 2 ) is the reduced mass. No spin-independent relativistic corrections are included. For the 
purpose of making some preliminary estimates of the energies of the lowest two states of the B c system, it is necessary 
to consider only the spin-spin part of the spin-dependent potential since these are S states. Thus, 

V SD -► V SS = -^^£ a (r)Si • Sa. (2) 
9mim2 

Our solutions to the Schrodinger equation are generated numerically with the given form for the central potential |l§| ] 
and the effects of the spin-dependent parts are added as a perturbation improvement. Thus, the lS-state hyperfine 
splitting is given by 

AE bb = r^-\R ls (0)\ a . (3) 

We address the question of the validity of first-order perturbation theory for the contact potential of Eq. (Q) in an 
Appendix. 

All of the potential parameters in this section are strictly flavor- independent. This includes any additive constants. 
The central potential parameters and the constituent mass parameters are fit to the low-lying energy levels of charmo- 
nium and the upsilon system. The strong coupling constant as is fit to the observed charmonium hyperfine splitting 
of 117 MeV. In Table I below we present results with 5 different central potentials. These include the Cornell potential 

El 

V(r) = Ar- n/r + C, (4) 
where A = 0.1756 GeV 2 , n = 0.52 and C = -0.8578GeV: Martin's power-law potential @, 

V(r) = -8.093 + 6.898 (r x lGeV) 01 , (5) 
(units of the potential are GeV), and the logarithmic potential 

V(r) = -0.6631 + 0.733 ln(r x lGeV), (6) 

(potential units are also GeV). Each of these forms was used by Eichten and Quigg. However, our values for the 
potential parameters are slightly different. We also have done calculations with the potential of Song and Lin |22j], 
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V(r) = Ar 1 / 2 - Br" 1 ' 2 , 
where A = 0.511 GeV 3 / 2 and S = 0.923 GeV 1 / 2 , and the Turin || potential, 

V(r) = -ar~ 3/4 + br 3/i + C, 



(7) 



(8) 



where a = 0.620 GeV 1/4 , b = 0.304 GeV 7/4 and C = -0.823 GeV. 

Our results in the first 3 columns of Table | are very similar to those presented in Table I of Eichten and Quigg [EJ . 
The results obtained with the Song-Lin and the Turin potentials in all cases fall between the extremes defined by the 
first 3 columns. Averaging over all 5 cases presented in Table H yields 



M Bc = 62511s 5 MeV, Mb* = 6328± 9 8 MeV, 



(9) 



as the predicted energies of the two lowest states. If we treat Eichten and Quigg's results in the same way, then we 
get 6258±i MeV for B c and 6333^° MeV for B*, in good agreement with our preliminary analysis. All of these 
results fall comfortably within the ranges found by Kwong and Rosner B, that is, 



6194 MeV < M Bc < 6292 MeV, 



(10a) 



6284 MeV < Mb* < 6357 MeV. 



(10b) 



III. RICHARDSON'S POTENTIAL AND THE SPIN-DEPENDENT POTENTIALS 



The starting point for the derivation p4| of Richardson's potential is the one-loop expression for the running 
coupling constant [§5j , namely, 
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(33-2n/)ln(|q 2 |/A a )' 



(11) 



where nj denotes the number of quark degrees of freedom accessible to the propagating gluon [g6| and A is the QCD 
scale. Richardson [Q] realized that one could tame the infrared singularity of Eq. ( pd| ) and produce a linear confining 
potential by making the replacement, 



|q 2 |/A 2 - |q 2 |/A 2 + l, 



(12) 



in the argument of the logarithm. Our 1991 calculation || was based on a modification of Richardson's potential 
suggested by Moxhay and Rosner pTjl , that is, the string constant is treated as a free parameter instead of being 
connected to the scale parameter A. Then the potential takes the form, 



V(r) = Ar 



(33-2n/)r 



where 
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ln(l + x 2 ) 
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To specify the spin-dependent potentials, we use a generalization of the Eichten-Feinberg fromalism |D| derived by 
Our notation is closely akin to that of Eichten and Quigg ||. Thus, 
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V 8 



D 
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2ml 



T' 



L-(Sj 
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L-(Si-S 2 ) . 
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S 
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T d , 



(15) 



where S12 denotes the tensor operator, that is, 4 [3Si • f S2 • f — Si • S2] and mi denotes the mass of the lighter 
(charmed) quark, if the two constituent masses are not equal. The quantities T are connected with the PTN potentials 
as follows: 
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2V 5 , T'=T a -4V 5 



(16a) 
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T h = - 
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T c = Va/12, T d = Va/B, 



(16b) 



where we have supplemented the PTN potentials with a long range spin-orbit contribution to include nonperturbative 
effects. Such a contribution is clearly indicated by the lattice calculations of Bali, Schilling and Wachter ]3| and is 
required to satisfy the consistency condition, 
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[E(r) + V 1 (r)-V 2 (r)]=0, 



(17) 



derived by Gromes |12] ] from the requirements of Lorentz covariance. 

The PTN potentials are expressed in terms of the strong coupling constant a§ and the renormalization scale /z. 
The function E(r) includes the leading correction to the central potential, 
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6ir 



(33-2n / )(ln/ir + 7£ ) + y-^- 



(18) 



where 7_e = 0.5772 • • •, is Euler's constant. This function is of course a part of the contribution to the short-distance 
behavior of the running coupling constant in Eq. ( |TT| ) , although its contribution to Richardson's potential is not made 
explicit. The first two spin-orbit potentials are given by 
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(19) 
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(33-2n / )(ln A ir + 7e ) + y - ^ - 9 (In y/rn^r + lE ) 



The tensor potential takes the form 
4a 



1 + a ^ 
6ir 
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(20) 
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The lengthy expression for the spin-spin potential 
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V 4 (r) = 
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reduces to an especially simple form for states with nonzero angular momentum, that is, 
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(22) 



(23) 



since one can take the derivatives present in Eq. (|22|) in this case, and the delta function there does not contribute. 
If the two constituent quarks have equal masses, then it is necessary to add a second-order contribution from the 
annihilation graphs, 



5F 4 = 8a|(l-ln2)<5 3 (r), 
to the expression of Eq. (|2|). The spin-orbit potential discovered by PTN is given by 

cA 

v 5(r) = — rln(mi/m 2 ). 



(24) 



(25) 



The equations describing the PTN potentials use the modified minimal subtraction scheme to define the strong 
coupling constant a§. This coupling constant can be related to that defined in a scheme developed by Gupta et al. 
Tol, which was used in the 1991 calculation S, that is, 
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to the second order in as- It is straightforward to verify that PTN's expressions for the potentials are the same as 
those of Gupta et al. to the second order in as by using Eq. ( p6|) to connect the two renormalization schemes. The 
signs of the terms in Eq. ( |2^ ) indicate that we should expect smaller coupling constants in the present calculation 
than in our 1991 calculation. 

In determining the matrix elements for the perturbation improvement of the central potential energies, we must 
exercise some care since the total spin, S = Si + S2, is not a good quantum number in the most general case. This 
is a consequence of the fact that its magnitude does not commute with the spin-dependent potential because, for 
example, 

[L-S!,5 2 ] = 2iL-(S 2 x Si). (27) 

We will deal with the more complicated case of the mass-mixing matrix in the next section. For the remainder of this 
section, we will consider only the diagonal matrix elements. This will suffice to determine the requisite expectation 
values for the case of equal masses, since then the total spin is a good quantum number. 

The diagonal matrix elements can all be evaluated in terms of the expectation values || of (L • S) since 

(L-Sr) = <L-S 2 > = ~<L.S>, (28) 
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and 



, s X _ 4(L2)(S2) - 6(L-S) - 12(L-S) 2 

{Sl2) (2/ + 3)(2Z-l) ' (29) 

where I is the orbital angular momentum quantum number, as we have verified [p8|-p0|. Hence, the diagonal matrix 
elements of the spin-dependent potentials may be expressed 



(V SD ) = <L • S) 



(T a ) , {T' a ) , (T b ) 



Amf 4771,2 vrt\vcii 



' [(S12) (T c ) + (Sr ■ S 2 ) (T d )} . (30) 
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Our results for the charmonium and upsilon energies obtained with Richardson's potential and the one-loop expres- 
sions for the spin-dependent potentials are presented in Tables || and [II (listed in column 4 as FUII98), where they 



are compared with the results of Eichten and Quigg || . The flavor-independent potential parameters and constituent 
masses used to obtain these results are 

A = 0.152 GeV 2 , A = 0.431 GeV, (31a) 



m b = 4.889 GeV, 



■m c 



1.476 GeV. 



(31b) 



For the upsilon system the value of the coupling constant ag — 0.30, 
fi = 1.95 GeV. For charmonium these values are = 0.486 and fi = 
universal QCD scale determined from 



A 



qcd = lie 



-67r/(33-2n f )ac 



and the value of the renormalization scale 
0.80 GeV. In both cases the value of the 



(32) 



is 0.190 GeV, a value consistent with other determinations |pjj| . Agreement between the calculated results and the 
measured values |3l]] in Tables || and III is extremely good, the average deviation being only 4.3 MeV for the upsilon 
system. Most of the 19.9 MeV deviation in the charmonium case arises from a 25 MeV discrepancy in the center-of- 
gravity of the IP state. Much of this difference could probably be removed by a proper account of the spin-independent 
relativistic corrections 11331 - 



Our results for the upsilon and charmonium lcptonic widths arc presented in Table IV, where they are compared 
with those calculated by Eichten and Quigg ||. Our decay rates were obtained with the formula fl33| ], 



M 2 (QQ) 



\R(0)\< 



16as 
3tt 



(33) 



where a is the fine structure constant, eg. denotes the quark charge, and M denotes the mass of the quark-antiquark 
state. Since the average values of the momentum transferred in the annihilation graphs underlying Eq. (|3^) should 
be much larger than that associated with the scattering processes of the spin-dependent potentials above, we should 
expect the appropriate value for as in Eq. (B3j) to be smaller than that used above for the fine structure splittings. 
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Thus we choose as = 0.18, a value obtained from heavy quarkonium decays ||g4|]. Our leptonic widths also agree with 
experiment very well. In order to carry out a fair comparison with the widths of Eichten and Quigg, it is necessary 
to correct their published results for radiative QCD corrections. Doing this as described above gives the EQ results 
shown in parentheses and moves their results much closer to the measured values. The differences between the values 



listed in columns 2 and 3 of Table IV gives a measure of the differences in wave functions calculated with Richardson's 
potential and the Buchmiiller-Tye potential. 

In order to have a good basis for comparison, we have carried out a Richardson's potential calculation with tree-level 
expressions for all the spin-dependent potentials. These may be readily obtained by omitting all terms of 0(a|) in 
the spin-dependent potentials above, and thus there is no reference to the renormalization scale. Then the expressions 
for the functions T are especially simple, that is, 

r. = 2*. = fH-£ (34a) 

n = |#, Tl = 0, T C =|§, T d = ^(r). (34b) 

The central potential parameters in this calculation are the same as those listed in Eq. (3l|). The constituent mass 
for the charmed quark is 1.476 GeV, the same as in Eq. (|3l|), and toj, = 4.884 GeV, slightly smaller. The coupling 
constant as = 0.339. Results of this calculation are also presented in Tables |0| and III (designated as FUI98 in column 



3). For both charmonium and the upsilon system, the average deviation is smaller for the one-loop calculation, but 
this does not tell the whole story. The fit to the fine structure splittings of charmonium with the one-loop calculation 
is much better than either of the tree-level calculations presented in Table O. This is our basis for the expectation 
that a one-loop calculation should give a more accurate rendering of the fine structure and the hyperfine structure 
of the B c system. Since the central potential parameters are the same, the leptonic widths for the FUI and FUII 
calculations are the same. 

One of the most dramatic differences in the one-loop and tree level results is in the prediction for upsilon hyperfine 



splittings in Table III. For example, the one- loop prediction for the IS state is 55 MeV, which is consistent with earlier 
predictions |8|,[l(| and the average of the tree-level predictions is 90 MeV. It is interesting to note that the one-loop 
level prediction is much closer to the lattice results (44-50 MeV) of Bali, Schilling and Wachter ||. Thus, measurement 
of the energies of the singlet S energies would serve to clarify the proper input for a good phenomenological calculation 
of the properties of heavy quark systems. Some experimental group should give these measurements a high priority. 

IV. MASS MIXING MATRIX FOR P STATES 

In determining the eigenvalues and eigenfunctions for the B c system, we have a choice of using basis functions from 
either the L-S coupling scheme or the j-j coupling scheme. We follow the lead of Eichten and Quigg || and choose 
the j-j basis. First we form the total angular momentum of the charmed quark, Ji = L + Si, and then we form the 
total angular momentum of the system, J = Ji + S2. For P states, the J = 2 states and the J = state are the same 
in either basis, that is, 

^2m(~~)= V>2m(ll), ^00(5 \) = Vbo(l 1), (35) 

and the determination of the eigenvalues of the spin-dependent potentials requires only the diagonal elements of Eq. 
(^0|). On the other hand, the J — 1 states are a combination, 

3 1 11 

^i m ^aii} lm {--)+a 2 ipi m {--), (36) 

and one must diagonalize the mass mixing matrix in order to determine the eigenvalues and eigenfunctions. In the 
basis described above, we have derived the following forms for the spin-dependent operators of Eq. (|l5|), that is, 




(L-Si) = I -1 I ■ (L • S 2 ) = I _±1 i ' j . 1 87a) 

and 




6 




2 2^2 

($,-,)= | a|l | | • (S, -S,i:: | £ _'l j . C-!7M 

To obtain the matrices for the operators L • S and L ■ (Si — Sa), one simply takes the appropriate linear combination 
of the matrices in Eqs. (37a). By inspection one can verify that our first 3 matrices in Eqs. ( |3^ ) are consistent with 
Eqs. (2.21) of Eichten and Quigg g. 

It is straightforward to take the equal-mass, or L-S, limit. Then one can find a unitary transformation that 
simultaneously diagonalizes the operators, L • S, Svz and Si • S2. The eigenvectors are given by 

AmCPl) = \fl^m(l \) + \j\^rn{\ \), (38a) 
VlmC^l) = t/f^mCf \) ~ J\^lm{\ \)- (38b) 



The first eigenvector of Eqs. (g8|) has the lower eigenvalue since it is determined by choosing the solution to the 
quadratic eigenvalue equation with the minus sign. Our choice of overall sign for the second eigenvector of Eqs. (|3^ ) 
is opposite that of Eichten and Quigg, but the same as that of Gershtein et al. M . 



V. RESULTS AND DISCUSSION 



Since the central potential parameters are strictly flavor independent and the constituent masses are not allowed to 
vary, the only decision that one has to make to address the B c system is to decide on the value of ag. The simplest 
choice is the average of the values in the charmonium and the upsilon system, that is, ag = 0.393. This choice requires 
a value of fx = 1.12 GeV, to preserve the value of 0.190 GeV for the QCD scale. Our results for the low-lying S, P 
and D states of the B c system are presented in column 6 of Table [v|, where they are compared with the predictions 
of Eichten and Quigg ||, Gershtein et al. Chen and Kuang |l6| ], the NRQCD lattice calculation |35| and our tree 
level calculation. The energy of each of these states lies below the B-D meson threshold at 7143 MeV. In determining 
the D-state splittings, which are much smaller than the P-state splittings, we followed the example of Eichten and 
Quigg and did not use the mass-mixing matrix approach. Several running coupling constant effects are noticeable in 
Table [v|, although our results are in general agreement with the earlier calculations. Our S-state hyperfine splittings 
are smaller than those of Eichten and Quigg, and our results for the P-state fine structure splittings are larger. We 
do not see the inverted order of the D states predicted by both Eichten and Quigg and Gershtein et al.. The order of 
the fine-structure levels |37],|38| is a consequence of the relative sizes of the the perturbative and the nonperturbative 
contributions to the spin-orbit potentials of Eq. (|l5|) . 

Comparison with the lattice NRQCD results ]3q ] is of special interest. Their result for the B c mass is based on 
a lattice calculation of the kinetic mass, 4.76 ± 0.02 and a determination of the scale factor by comparison with 
charmonium and the upsilon system, which gives a -1 = 1.32 ± 0.04 GeV. Thus, we list two errors for the lattice 
calculation of the B c mass in Table [v| The smaller error is a consequence of the error in the kinetic mass calculation 
and the larger error is a consequence of the uncertainty in the overall scale. In determining the other lattice errors, 
we simply list the largest possible source of error without including the error due to the overall scale. The lattice 
result for the IS hyperfine splitting is 41 ± 3 MeV, which is independent of the larger errors discussed above, and 
somewhat smaller than our 55 MeV result. The lattice result for the IP 2 — IP splitting is 59 ±5 MeV, in reasonable 
agreement with our result of 71 MeV. Although the number quoted in Ref. J35| for the B c mass is very close to our 
result, its importance as a confirmation of our work is undermined by the large error in the overall scale. Another 
important lattice result is now available from the UKQCD collaboration |36|. Their result for the ground state energy 
is 6386±9±98±15 MeV, which is consistent with our result because of the large error bars. It will be interesting to see 
whether further refinements of the lattice calculations will support our results, or offer the experimenters alternative 
predictions. 

In order to get some idea of an error estimate for our predictions, we have calculated the ground state mass of the 
B c system as a function of a§ in a range bounded by the values determined in charmonium and the upsilon system. 
Our results are shown in Fig. [lj where they are compared with the result of Table |, the Eichten-Quigg quote and 
the NRQCD lattice result. The error shown for the lattice result simply ignores the large error associated with the 
overall scale. The horizontal lines there describe the limits determined by Kwong and Rosner [Q. Using the largest 
and smallest values of ag in Fig. [l] to determine the errors, we have 
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M Bc = 6286±£ 5 MeV M B , = 6341±| MeV, (39) 

as our predicted value for the IS energies. It is interesting to note how close our results are to the earlier predictions 
of Godfrey and Isgur JIJ] (6270 and 6340 MeV) and the predictions of Baker, Ball and Zachariasen (6287 and 
6372 MeV), although our prediction for the hyperfinc splitting is smaller than either. Clearly, the precision of the 
experiments |lj requires a very substantial improvement to be sensitive to the energy differences between the various 
calculations listed in Table [v| and Fig. [p. 

The pseudoscalar decay constant is given by the Van Royen-Weisskopf formula modified for color |39||, that is, 



Ik = ,40) 



and we find that 

f Bc = 517 MeV, (41) 
in excellent agreement with Eichten and Quigg's result, f Bc = 500 MeV, and in reasonable agreement with the lattice 



result |35|, f Bc = 440 ± 20 MeV. 



The empirical result obtained by Collins et al. J40| for potential model wave functions at the origin, that is, 

I^WI^I^/^o)! 1 - 3 !^^)! - 7 , (42) 

provides another touchstone for our numerical work. Using this relationship and input from our charmonium and 
upsilon calculations, we get |i?e e (0)| 2 = 1.81 GeV 3 , about 3% higher than our numerical result. Using their empirical 
relationship for the ground state hyperfine splittings, that is, 

M B , - M Bc ~ 0.7 (M J/4 , - M„ c )°' 6S {Mr ~ M Vb ) ' 35 , (43) 

yields a splitting of 63 MeV, about 14% larger than our result listed in Table V. Both of these results are reasonable 
in view of the spread of the results Collins et al. obtained with different forms for the central potential. 

Diagonalizing the P-state mixing matrix, we obtain the following combinations for the two lowest J = 1 P states, 

^i m (l + ) = 0.118 i/j lm {- i) + 0.993 Vim(^ \), (44a) 

V-im(l + ') = 0.993 Vi m (| i) - 0.118Vim(i i), (44b) 

which is very close to the j-j coupling limit, that one would expect to be valid in the heavy-quark limit. Using the 
inverse of Eqs. (|§8|), we can determine the probability of observing spin 1 in the lowest 1 + state, P x + (5 = 1) = 0.773. 
Our result is consistent with the lattice calculation |35fl , where the mixing angle in the L-S basis was found to be 
close to that of the j-j limit. From this mixing angle Hd = 33.4° ± 1.5°) we obtain P 1+ (S = 1) = 0.697 ± 0.020, in 
reasonable agreement with our result. Our results for the mixing angles of the 1 + and 1 + states is very different from 
that of Eichten and Quigg, whose results were much closer to the L-S limit. Below we show that this difference has 
important implications for the spectrum of photons emitted in electric dipole transitions. 
The electric dipole rate for the emission of a photon JsJ of energy k is given by 

T E i(i - / + 7) = --^fc 3 (2J / + 1)1 (/|r|») | 2 %, (45) 

where the statistical factor 5j/ = 1 for transitions between triplet S and triplet P states H] , and 5,/ = 3 for transitions 
between spin-singlet states. The mean charge in Eq. |l|) is 

m 2 ei-mie 2 , . 

( e Q> = . (46) 

mi + TO2 

where e% denotes the charge of the charmed quark (in units of the proton's charge) and e 2 denotes the charge of the 



bottom antiquark. Our results for the IP — > 15 and the 25 — * IP transition rates are shown in Table VI, where they 
are compared with Eichten and Quigg's results. Our rates for transitions involving the IP 2 and IP states are rather 
close to their counterparts calculated by Eichten and Quigg. However for the J = 1 states important differences arise. 
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In particular, we predict 4 nonzero transition probabilities instead of 2. Thus, each of our simulated photon spectra 
presented in Fig. g and Fig.f3 has six lines instead of 4. These additional lines are a consequence of the fact that our 
mixture of states is in Eq. ([44) is not close to the L-S limit, and thus both 1 + states have a substantial admixture of 
both triplet and singlet components. 

The magnetic dipole transition rate between S states is given by 

r M i(i - / + 7 ) = ^L 9 fc 3 (2J/ + l)|</bo(fcr/2)N)| 2 , (47) 



where the mean magnetic dipole moment is 



m 2 e 1 -m 1 e 2 

Mmag = ■ (48) 

4771x7712 



Our results for the magnetic dipole transition rates are presented in Table |VII| , where they are compared with those 
of Eichten and Quigg || and Gershtein et al. M . Most of the differences are a consequence of different energies for 
the hyperfine splittings, since the results for the matrix elements are rather close. 



The photon energies and transition rates in Tables VI and VII suggest at least two good strategies for experimental 
searches for the IP and IS states. One could detect one or both of the high energy photons (457, 436 MeV) associated 
with the decay of the lowest two 1 + states to the ground state and then search for a leptonic decay of the ground 
state. An alternative would be to look for some of the high-energy photons (417, 406, 384, 350 MeV) in the IP — > l 3 5i 
transitions and then seek a coincidence with the 55 MeV photon associated with the decay to the ground state. 

Since the charmed quark and the bottom antiquark cannot annihilate into gluons, the only additional complication 
that arises in the decays of these low-lying states is the tt tt channel. If we take the rates for the tt tt decays from 
Eichten and Quigg, then we can work out complet e dec ay schemes and branching ratios for the 2S and IP states as 



well as the l 3 5i state. These are shown in Table VIII . Our table of decay rates and branching ratios differs from 
that of Eichten and Quigg in two important respects. More photon channels are available to the 1 + states and the 
2S states. Our decay width for the 1 3 S\ is 59 eV, more than a factor of 2 smaller than theirs. 

Our final calculation is that of the lifetime of the ground state of the B c system. We follow the approach used 
by several researchers |f42|-[4"4|| where the decay of the B c meson is taken to be the sum of 3 distinct contributions, 
namely weak decay of the b antiquark while the c quark behaves as a spectator, weak decay of the c quark while the b 
antiquark behaves as a spectator and an annihilation of the b antiquark and the c quark into an intermediate vector 
boson that subsequently decays into a lepton- neutrino pair or a quark- antiquark pair. Thus the total decay rate is 
the sum, 

T(B C -> X) = T(b -> X) + r(c -> X) + T(annih). (49) 
If one neglects quark binding effects, then the first two terms are given by 

r( ^ x) = 52i|«, r( _ x) = W^, (50) 

1927T' 3 1927r J 

where the subscripted quantities V denote the appropriate Cabibbo-Kobayashi-Maskawa matrix elements EH] and Gf 
denotes the Fermi coupling constant. Using the constituent masses listed in Eq. (|l|), we obtain 

r(6 -+X) = (8.73 ± 1.34) x 10" 10 MeV , T(c X) = (7.59 ± 0.02) x 10~ 10 MeV. (51) 

The annihilation width is given by 

T(annih) = %V bc \ 2 f B M Bc £ m? (l - H-) C h (52) 

i ^ B c / 

where m, denotes the mass of the heavier Fermion in the given decay channel. The most important channels in the 
sum are the rv and the cs channels. For the former Cj = 1, and for the latter Ci = 3\V CS \ 2 . Using Eq. (Ell) for /s c 



and Eq. (39) for Mg c , we have 

T(annih) = (1.13 ± 0.17) x 10~ 10 MeV. (53) 
Adding these 3 widths yields a lifetime 

t Bc = 0.38 ±0.03 ps, (54) 
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in good agreement with the measured CDF result m. Our result is also in reasonable agreement with the recent 
calculation of El-Hady, Lodhi and Vary [jjij , who obtain tb c — 0.46 ps, since the spectator widths of Eqs. j50| ) 
are very sensitive to small differences in the constituent masses. Although different authors may wish to interprete 
relatively small differences between T(b — > X) and T(c — > X) as the domination of one process over the other, we feel 
that the safest characterization of our results is that the two spectator processes are almost equally important and 
that the annihilation processes are less important, consistent with the earlier conclusion of Gershtein et al. |4^| . 
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APPENDIX A: HYPERFINE SPLITTING AND CONTACT POTENTIALS 



Determining the energy shifts from the hyperfine splitting can be problematic because of the delta function that is 



often present in the spin-spin potential. As Lucha, Schoberl and Gromes point out in their review |24|, the energy of 
the singlet state is not actually bounded from below, in marked contrast to the first-order perturbation theory result 
ofEq. (g). Although the use of Eq. (g) in the literature is fairly common, seldom does one see any discussion of its 
validity. It is straightforward to create a context for addressing this question, by considering a more general form for 
the spin-spin potential, which allows for a finite range, that is, 

Vss = Q 32a ° "Tla Si-Sa- (Al) 

The advantage of such a form is that one can calculate the singlet and triplet S-state eigenvalues of the Hamiltonian 
of Eq. (0) exactly and examine the limit as the range parameter b becomes smaller and smaller. We have done 
such a calculation with the logarithmic potential of Eq. (g) and the parameters listed in Table I. Our results for the 
exact singlet and triplet as a function of the range b are shown in Fig. 0, where they are compared with the results 



of a first-order perturbation calculation of the singlet and triplet energies produced by the potential of Eq. (Al) 
The contact potential results of Eq. (g) are presented as two horizontal lines. It is gratifying to see the first-order 



perturbation result from Eq. ( |Al| ) above approach the contact potential result in the limit 6 — > 0. Extrapolating the 
exact result for the triplet energy to the 6 = limit gives 6330 MeV, about 3 MeV lower than the result listed in 
column 4 of Table 1, which was obtained with the contact potential. Thus, it is clear that first-order perturbation 
theory and the contact potential give a good account of the IS triplet energy. 

Fig. [| gives a clear signal of the instability of the singlet energies as b — > 0. However the difficulty begins to 
appear only as b decreases below 0.3 GeV -1 . The short range required for the effects of the instability to manifest 
itself provides a means of resolving this dilemma. As Lucha, Schoberl and Gromes point out, the contact potential 
expression is not really valid for such short ranges since one must take the nonrelativistic limit in order to obtain 
it. Such a limit requires that nonlocal effects associated with the normalization of Dirac wave functions be ignored. 
Thus it might be reasonable to choose b — 0.36 GeV, the geometric mean of the Compton wavelengths of the two 
constituent quarks. Such a choice would lead to a value of the singlet energy very close to that found in column 4 of 
Table 1. 
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FIG. 1. Ground state energy of the B c system as a function of the running coupling constant. The PH98 result is based 
the calculations presented in Table || 

FIG. 2. Simulated photon spectrum for IP — > IS* transitions. The probability of populating each of the initial states 
assumed to be equal. 

FIG. 3. Simulated photon spectrum for 2S — > IP transitions. 
FIG. 4. Hyperfine splittings of the IS state for a finite range spin-spin potential. 
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TABLE I. Ground state energies of the heavy-quark systems (MeV). 



State 


Cornell 


Martin 


logarithm 


Song-Lin 


Turin 


as 


0.313 


0.437 


0.372 


0.396 


0.373 


m b (GeV) 


5.232 


5.174 


4.905 


5.199 


5.171 


m c (GeV) 
l 3 5i (cc) 


1.840 


1.800 


1.500 


1.820 


1.790 


3097 


3097 


3097 


3097 


3097 


1% 


2980 


2980 


2980 


2980 


2980 


AE 1S 


117 


117 


117 


117 


117 


l 3 ^! (66) 


9461 


9461 


9460 


9460 


9460 


1% 


9316 


9397 


9395 


9380 


9365 


A£i S 


145 


64 


65 


80 


95 


l 3 Si(c6) 


6337 


6319 


6333 


6324 


6327 


1% 


6246 


6247 


6266 


6247 


6247 




91 


72 


67 


77 


80 



TABLE II. Charmonium energies (MeV). 



State 


EQ94 


FUI98 


FUII98 


EXPT 


1% 


3097 


3098 


3098 


3097 ± 0.04 


l 1 ^ 


2980 


2981 


2980 


2980 ± 2 


2 3 Si 


3686 


3692 


3693 


3686 ± 0.1 


2 1 S 


3608 


3617 


3615 




1 3 P 2 


3507 


3515 


3530 


3556 ± 0.1 


l 3 Pi 


3486 


3492 


3482 


3511 ± 0.1 


l 3 Po 


3436 


3443 


3391 


3415 ± 1.0 


l^i 


3493 


3499 


3501 


3526 ± 0.1 


VP 


25.5 


22.7 


19.9 





TABLE III. Upsilon energies (MeV). 



State 


EQ94 


FUI98 


FUII98 


EXPT 


l 3 Si 


9464 


9461 


9461 


9460 ± 0.2 


1^0 


9377 


9368 


9406 




2 3 Si 


10007 


10022 


10027 


10023 ± 0.3 


2 1 S*o 


9963 


9978 


10001 




3^1 


10339 


10358 


10364 


10353 ± 0.5 


3^0 


10298 


10325 


10344 




1 3 P 2 


9886 


9902 


9910 


9913 ± 0.6 


l 3 Pi 


9864 


9881 


9891 


9892 ± 0.7 


l 3 Po 


9834 


9852 


9863 


9860 ± 1 


l^i 


9873 


9889 


9899 




2 3 P 2 


10242 


10261 


10269 


10269 ± 0.4 


2 3 Pi 


10224 


10244 


10255 


10255 ± 0.5 


2 3 P 


10199 


10222 


10234 


10232 ± 0.6 


2 x Pi 


10231 


10251 


10261 




VP 


24.9 


8.3 


4.3 
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TABLE IV. Leptonic widths (KeV). 



State EQ94 FU98 EXPT 

T(15) 1.71 (1.19) L34 1.32 ± 0.05 

T(2S) 0.76 (0.53) 0.57 0.52 ± 0.03 

T(35) 0.55 (0.38) 0.40 0.48 ± 0.06 

tp(lS) 8.00 (5.55) 5.81 5.26 ± 0.37 

ip{2S) 3.67 (2.55) 2.61 2.14 ± 0.21 



TABLE V. Energies of the B c system (MeV). 



State 


EQ94 


GKLT95 


CK92 


FUI98 


FUII98 


LAT96 


l 3 Si 


6337 


6317 


6355 


6341 


6341 


6321 ± 30 


1% 


6264 


6253 


6310 


6267 


6286 


6280 ± 30 ± 190 


2 3 Si 


6899 


6902 


6917 


6911 


6914 


6990 ± 80 


2 1 S 


6856 


6867 


6890 


6869 


6882 


6960 ± 80 


1P 2 


6747 


6743 


6773 


6761 


6772 


6783 ± 30 


1P1+' 


6736 


6729 




6750 


6760 


6765 ± 30 


1P1+ 


6730 


6717 




6742 


6737 


6743 ± 30 


IPo 


6700 


6683 


6728 


6713 


6701 


6727 ± 30 


1 3 P>3 


7005 


7007 




7022 


7032 




1 3 P»2 


7012 






7025 


7028 




l 3 Dl 


7012 






7024 


7019 




1^2 


7009 


7008 




7023 


7028 





TABLE VI. Electric dipole matrix elements and transition rates. 



Transition 


Photon energy 
(MeV) 


</|r|i> (GeV)- 1 
GKLT95 EQ94 


FU98 


r(KeV) 
EQ94 


FU98 


1P 2 -> l 3 Si 


417 






113 


126 


IP 1+' -> l 3 Si 


406 






0.1 


26.2 


IP 1+ -> l 3 Si 


384 






99.5 


75.8 


IPo l 3 Si 


350 


1.568 1.714 


1.683 


79.2 


74.2 


IP 1+' -> 1% 


457 






56.4 


128 


IP 1+ l^o 


436 






0.0 


32.5 


2 3 Si 1P 2 


141 






17.7 


14.5 


2 3 Si -> IP 1+' 


152 






0.0 


2.5 


2 3 Si -> IP 1+ 


175 






14.5 


13.3 


2 3 Si IPo 


210 


-2.019 -2.247 


-2.253 


7.8 


9.6 


2 1 So -> IP 1+' 


121 






5.2 


13.1 


2 1 5 -> IP 1+ 


143 






0.0 


6.4 



TABLE VII. Magnetic dipole matrix elements and transition rates. 



Transition 


Photon energy 
(MeV) 


</|jb(fcr/2)|t) 
EQ94 


FU98 


GKLT95 


T(KeV) 
EQ94 


FU98 


2 3 5i -> 2% 


32 


0.9990 


0.9995 


0.010 


0.029 


0.012 


2 3 S'i -> l^o 


599 


0.0395 


0.0399 


0.098 


0.123 


0.122 


2^0 l 3 S*i 


520 


0.0265 


0.0305 


0.096 


0.093 


0.139 


l 3 Si - l^o 


55 


9.9993 


0.9996 


0.060 


0.135 


0.059 
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TABLE VIII. Decays and branching ratios of the B c system. 



State 


Total width (KcV) 


Decay mode 


Branching ratio [ 


rSi 


0.059 


1 So ± 7 


100 


m 


126 


l 3 Si + 7 


100 


ipi+' 


154 


l 3 Si + 7 


17 






l'So + 7 


83 


IP i+ 


108 


l 3 Si ± 7 


70 






l'So + 7 


30 


lPo 


74.2 


iis-i ± 7 


100 


2 3 Si 


89.9 ±7 


l 3 Sl + 7T7T 


56 ±8 






1P 2 ±7 


16 ± 1 






1P1 + ' +7 


3 ±0.2 






1P1+ ±7 


15 ± 1 






lPo + 7 


11 ±1 


2 1 S 


69.5 ±7 


1 So ± VT7T 


72 ±10 






1P1+' +7 


19 ±2 






1P1+ +7 


9± 1 



(percent) 



14 



o 

CO 



E— 1 
< 
E— 1 

GO 



J I I I I I I I I I I I I I I I I I I I I I I I I I L 



O 

CO 



o 



^ CD 



O ?H 

o a; 

CO 

o 

°^ 

CO ft 



o 

CO 
CO 



o 

CO 



o 

cv 

o o o o o o o o 00 

^ C\2 O CO CO ^ C\2 

T 1 T 1 T 1 



o 

C\2 
C\2 



E— 1 

<J O 

E— 1 — O 

GO C\2 



00 o 
C\2 —I CO 



> 

CD 



O ?H 

CO CD 

t — i d 

CD 

o 
■+J 
o 

o 



o 

C\2 



J L I I I I I I I I I I I I I I I I § 

LO O LO O ^ 




(A90) sai§jaua si 



